%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Square well potential.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function QA_Movimento_SquareWell
	figure(1);					clf;

	uA=1.66e-27;				m =87*uA;									%atomic mass unit & mass of strontium
	L = 1e-6;					x =[0:.005:1]*L;							%size of potential
		
	for nn=1:6
		E =nn^2;
		ps=(sqrt(2/L)*sin(nn*pi*x/L)*1e-3).^2 + E;
		subplot(121),fill(x/L,ps,'b','LineWidth',1.4)
		set(gca,'xLim',[0 1],'FontSize',16,'LineWidth',1.4);%'yLim',[0 E],
		xlabel('\fontsize{16}{\itx / L}   (\mum)');
		ylabel('\fontsize{16}{\itE} , {\it\psi}   (\mum)');
		hold on;
	end
	subplot(121),plot(x*1e6,ones(size(ps)),'r','LineWidth',1.4)
end